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DIRECTIONALITY  OF  CONTINUUM  GAMMA  RAYS 


FROM  SOLAR  FLARES 
TAEIL  BAI 

Center  for  Space  Science  and  Astrophysics,  Stanford  University 

ABSTRACT 

Using  the  hard  X-ray  spectrometer  total  counts  as  calibration.  I 
investigate  how  the  brightness  of  y-rays  above  300  keV  changes  as  a 
function  of  the  heliocentric  angle  of  a  flare.  The  normalized  y-ray 
brightness,  on  the  average,  increases  with  the  heliocentric  angle;  a  flare  at 
the  limb  is  estimated  to  be  13  times  brighter  in  y-rays  than  a  similar  flare 
near  the  central  meridian.  Both  pancake-like  electron  distributions  and 
downward  beam  distributions  can  be  adjusted  to  produce  the  deduced  limb 
brightening. 


I.  INTRODUCTION 

Radiation  from  energetic  electrons  is  directional,  with  more 
radiation  being  emitted  in  the  forward  direction.  The  more  energetic  the 
electron,  the  more  directional  the  radiation  becomes.  In  solar  flares  some  of 
the  hard  X-ray  photons  emitted  toward  the  photosphere  are  backscattered 
via  Compton  scattering  and  the  rest  are  consumed  in  photoionization 
(Langer  and  Petrosian  1977;  Bai  and  Ramaty  1978).  Backscattered  photons 
tend  to  diminish  the  directionality  of  observed  hard  X-rays;  hence,  the 
directionality  of  hard  X-rays  may  be  detectable  only  at  high  energies  >100 
keV,  depending  on  the  distribution  of  electrons  (Bai  and  Ramaty  1978; 
Leach  and  Petrosian  1983).  Therefore,  the  anisotropy  of  hard  X-rays  and  y- 
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rays  from  solar  flares  can  yield  information  on  the  angular  distribution  of 
the  momenta  of  energetic  electrons. 

There  are  several  ways  to  study  the  anisotropy  of  solar  hard  X-  and  y- 
rays.  One  is  to  study  the  distribution  on  the  disk  of  the  Sun  of  flares  that  are 
detected  with  a  certain  detector.  Suppose  that  more  radiation  is  emitted  in 
the  directions  parallel  to  the  surface  of  the  Sun  than  in  the  vertical 
dneouun.  Then,  a  flare  located  near  the  limb  will  look  brighter  than  the 
same  flare  located  near  the  disk  center.  This  effect  is  called  "limb 
brightening."  Limb  brightening  will  in  turn  push  above  the  detection 
threshold  many  limb  flares  which  would  not  be  detected  were  they  located 
near  the  disk  center.  Distributions  of  flares  detected  with  hard  X-ray 
detectors  have  been  studied  (Kane  1974;  Datlowe  et  al.  1977).  But  because  of 
the  effect  of  Compton  backscatter,  the  limb  brightening  is  negligible  at  the 
photon  energies  in  their  studies,  even  if  the  electron  distribution  is 
anisotropic  (Bai  and  Ramaty  1978;  Langer  and  Petrosian  1977;  Leach  and 
Petrosian  1983).  Hence,  one  could  not  draw  any  firm  conclusions  from  their 
studies.  However,  at  higher  energies  the  anisotropy  is  detectable,  if  the 
electron  distribution  shows  a  reasonable  degree  of  anisotropy.  Vestrand  et 
al.  (1987)  studied  the  distribution  of  flares  detected  with  the  gamma  ray 
spectrometer  (GRS)  aboard  the  Solar  Maximum  Mission  (SMM)  whose 
threshold  energy  is  300  keV.  They  found  that  more  flares  are  detected  near 
the  limb  than  are  expected  from  an  isotropic  radiation  case.  The  excess  of 
limb  flares  is  at  the  2.5  a  level.  Such  limb  brightening  can  be  caused  either 
by  electrons  streaming  down  toward  the  photosphere  or  by  electrons  with 
pancake-like  momentum  distributions  (Dermer  and  Ramaty  1986). 

The  second  way  of  studying  the  momentum  distributions  of  energetic 
electrons  is  to  investigate  photon  spectra.  Because  the  anisotropy  of 
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radiation  increases  with  increasing  photon  energy  for  a  given  electron 
momentum  distribution,  the  photon  spectrum  of  a  given  flare  will  vary 
depending  on  its  location  on  the  Sun.  By  studying  the  y-ray  continuum 
spectrum  of  the  GRS  events,  Vestrand  et  al.  (1987)  found  that  the  y-ray 
spectrum  on  the  average  becomes  flatter  as  the  location  of  flares  is  closer  to 
the  limb.  This  effect  is  significant  at  a  3.4  a  level.  The  observed  spectral 
behavior  is  also  consistent  with  either  downward  beam  distributions  or 


pancake-like  momentum  distributions. 

The  third  way  of  studying  the  momentum  distributions  of  energetic 
electrons  is  to  observe  the  same  flares  at  different  viewing  angles.  If  we 
have  two  identical  hard  X-ray  detectors  observing  the  Sun  at  widely 
separate  angles,  the  task  of  unraveling  the  radiation  pattern  will  be  rather 
simple.  But  stereoscopic  observations  analyzed  by  Kane  et  al.  (1988)  were 
made  with  two  different  detectors  with  different  efficiencies  and  energy 
channels.  Therefore,  a  large  number  of  flares  are  needed  to  compensate  for 
such  deficiency.  Kane  et  al.  (1988)  analyzed  39  flares  commonly  observed 
without  limb  occultation  by  the  hard  X-ray  detectors  aboard  the  Pioneer 
Venus  Orbiter  (PVO)  and  the  third  International  Sun-Earth  Explorer  (ISEE 
3;  now  ICE).  According  to  Kane  et  al.  (1988),  the  observed  anisotropy  is 
much  smaller  than  what  is  expected  from  model  calculations  made  by 
Petrosian  (1985)  and  Dermer  and  Ramaty  (1986)  to  fit  the  GRS  flare 
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In  order  to  help  resolve  this  conflict,  I  have  made  the  following  study. 
Because  of  Compton  backscatter,  the  directionality  of  30  keV  photons  is 
negligible,  even  if  the  electron  distribution  is  highly  anisotropic  (Bai  and 
Ramaty  1978;  Leach  and  Petrosian  1983).  Since  the  threshold  energy  of  the 
hard  X-ray  burst  spectrometer  (HXRBS)  is  about  30  keV,  we  can  use  the 
HXRBS  total  counts  as  a  calibrator  in  evaluating  the  directivity  of  y-ray 
continuum  emission. 

Figure  1  shows  the  correlation  between  HXRBS  total  counts,  X,  and 
GRS  fluences,  Y.  for  flares  observed  with  GRS  aboard  SMM  between  1980 
February  and  1985  December.  GRS  fluences  are  time  integrated  flux  of  y- 
rays  above  300  keV.  GRS  fluences  are  evaluated  from  the  formula 

y=A  0.3-1*-') I (8-1)  (1) 

where  A  is  the  differential  y-ray  fluenee  at  1  MeV  and  5  is  the  spectral  index 
given  in  Table  1  of  Vestrand  et  al.  HXRBS  total  counts  are  taken  from 
Dennis  et  al.  (1985).  Solid  dots  in  Figure  1  indicate  disk  flares,  and  "  +  " 
marks  indicate  limb  flares.  To  be  consistent  with  Vestrand  et  al.,  flares 
with  heliocentric  angle  greater  than  64  degrees  are  defined  as  limb  flares. 

In  order  to  determine  the  limb  brightening  of  y-rays  above  300  keV, 
one  can  make  least  square  fits  to  the  data  points  in  Figure  1.  The  GRS  flares 
are  selected  by  virtue  of  emitting  y-rays  above  the  GRS  threshold.  Because 
HXRBS  total  counts  and  GRS  fluences  have  a  positive  correlation,  many 
flares  with  small  HXRBS  total  counts  did  not  emit  detectable  y-rays.  To 
remove  the  bias  of  the  data  set  due  to  the  GRS  threshold  effect,  I  first 
selected  flares  with  HXRBS  total  counts  greater  than  6x10^  in  making  the 
fits.  There  are  60  flares  meeting  this  requirement.  Second,  even  among 
flares  with  HXRBS  total  counts  greater  than  6x10^  some  were  not  detected 
by  the  GRS.  To  compensate  for  this,  I  weighted  the  flares  with  GRS  fluences 
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below  10  photons  cm~2  by  the  inverse  of  the  detectability  function.  The  y-ray 
size  distribution  function,  N(Y),  which  shows  the  number  of  flares  above  a 
given  y-ray  fluence  Y,  is  well  fitted  by  a  power  law  down  to  Y=10  photons 
cm~-  (e.g.,  Chupp  1984).  The  deviation  from  the  power  law  below  10 
photons  cm-2  is  assumed  to  be  caused  by  the  detectability. 

After  these  two  procedures  I  determined  parameters  a  and  P  which 
minimize  the  quantity 

n 

S=^(yi  -a-p*f-)2  (2) 

i=l 

where  x  £-  =  logw  (Xt!  1(P)  and  y  i=log  10Y i,  X{  is  the  HXRBS  total  count,  and 
Yj  is  the  GRS  fluence  of  the  ith  flare.  Parameters  a  and  P  are  determined 
for  disk  flares  and  limb  flares  separately.  The  linear  regression  line  thus 
obtained  for  disk  flares  is 

yd=1.04x,  (3) 

and  it  is  shown  by  the  solid  line.  The  linear  regression  line  for  limb  flares  is 

y,-0.74+0.92x,  (4) 

and  is  shown  by  the  dashed  line.  Notice  here  that  the  latter  is  located  higher 
than  the  former.  For  x=l,  Ay=y/ -yd=0.62\  this  means  that  at  X=106  counts 

limb  flares  are  brighter  in  y  rays  by  a  factor  f=10°-62=4.2.  Similarly,  at 
,Y=  If)7  counts  limb  flares  are  brighter  in  y  rays  by  a  factor  of  3.2.  The 
change  of  the  brightening  factor  with  X  is  due  to  the  difference  in  slopes  of 
equations  f3)  and  (4).  However,  this  difference  is  not  statistically  significant. 
Therefore,  we  can  say  that  limb  flares  are  brighter  in  y-rays  than  disk 
flares  by  about  3.5  on  the  average. 

One  can  make  least-squares  fits  by  minimizing  the  quantity 
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(5) 


S*=i  (Xi-C-Dy,)2 

i=l 

The  regression  line  obtained  with  equation  (2)  can  be  used  to  estimate  the 
expected  y-ray  fluence  for  a  given  value  of  HXRBS  total  counts.  On  the  other 
hand,  the  regression  line  obtained  with  equation  (5)  can  be  used  to  estimate 
the  expected  HXRBS  total  counts  for  a  given  value  of  y-ray  fluence. 
Therefore,  equation  (2)  is  physically  more  meaningful  in  determining  the  y- 
ray  brightening  factor.  However,  equation  (5)  is  more  advantageous  in  that 
we  do  not  have  to  worry  about  the  GRS  threshold  effect.  The  linear 
regression  line  thus  obtained  for  disk  flares  is 

yd=0.20+1.18x,  (6) 

and  that  for  limb  flares  is 

yi=0.06+1.61x.  (7) 

From  equations  (6)  and  (7)  we  estimate  the  limb  brightening  factor  of  y  rays 
to  be  2.0  for  x=l  Qf=106  counts)  and  5.2  forX=107  counts.  These  estimates  are 
in  rough  agreement  with  those  obtained  from  equations  (3)  and  (4). 

I  have  done  a  similar  analysis  by  using  HXRBS  peak  rates  instead  of 
HXRBS  total  counts.  The  results  are  very  similar. 

Not  only  are  the  GRS  fluences  and  HXRBS  total  counts  positively 
correlated  but  also  the  coefficients  of  the  regression  lines  given  by  equations 
(2)  and  03)  are  almost  equal  to  unity.  Therefore,  we  can  with  some 
confidence  normalize  the  GRS  fluences  by  dividing  by  their  HXRBS  total 
counts.  I  define  the  quantity  k  as 

k  =  logi0(105Y/X)  (8) 

where  Y  is  the  GRS  fluence  above  300  keV  and  X  is  the  HXRBS  total  count. 
This  quantity  is  plotted  against  the  heliocentric  angle  0  in  Figure  2.  The 
straight  line  is  the  best  fit  to  the  data,  which  is  expressed  by 
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k  =  0.27  +  0.0147(6-61).  (9) 

The  correlation  coefficient  r  between  6  and  k  is  0.39  for  60  data  points.  (Here 
I  used  the  same  criterion  for  the  selection  of  flares  and  the  same  weighting 
as  in  obtaining  the  regression  lines  shown  in  Fig.  1.)  The  chance  for  60 
pairs  of  uncorrelated  values  to  have  r>0.39  is  only  0.2%.  The  uncertainty  in 
the  estimate  of  the  coefficient  0.0147  is  0.0041;  therefore,  the  hypothesis  that 
y-ray  emission  above  300  keV  is  not  directional  is  ruled  out  at  the  3.6  o  level 
0.0147/0.0041=3.6).  The  y-ray  brightening  factor  increases  on  the  average  as 
flares  move  out  from  the  disk  center  to  the  limb.  A  flare  at  the  limb  is 
expected  to  be  brighter  in  y  rays  on  the  average  by  a  factor  f  =  10°-<)147(90-15)  = 
12.7  than  the  same  flare  at  the  central  meridian  (1-0°,  b=15°).  (The 
heliocentric  angle  is  given  by  the  expression  9=  arccos  (cos  l  cos  b),  where  l 
is  the  longitude  and  b  is  the  latitude.  Because  flares  seldom  occur  near  the 
equator ,  few  flares  have  lvalues  less  than  15  degrees.) 

The  center-to-limb  variation  function  f  can  be  obtained  from  equation 
(9): 

f(6)  =  10  0.0147(0-51)'  (10) 

where  f  is  normalized  to  1  at  6  =  51°,  the  mean  heliocentric  angle  of  the 
flares  used  in  the  analysis.  If  we  remove  the  center-to-limb  variation  by 
dividing  GRS  fluences  by  f(6),  the  correlation  coefficient  between  GRS 
fluences  and  HXRBS  total  counts  increases  from  0.49  to  0.54.  The 
correlation  diagram  between  the  GRS  fluences  divided  by  fI9)  and  HXRBS 
total  counts  is  shown  in  Figure  4.  We  find  that  Figure  4  shows  much  less 
scatter  than  Figure  1.  The  scatter  in  Figure  4  is  mostly  due  to  variability  of 
individual  flares.  The  electron  spectrum  varies  from  flare  to  flare,  and  the 
anisotropy  of  energetic  electrons  may  vary  from  flare  to  flare. 
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Dermer  and  Ramaty  (1986)  calculated  the  directionality  of 
bremsstrahlung  of  energetic  electrons  for  various  electron  distributions 
described  by  the  following  two  distribution  functions.  The  first  represents 


downward  beam  distributions: 

M /cos  a)  = 


r  exp(r  cos  a) 


(10) 


2n[exp(r)-  exp(-  r)J 
whei-e  a  is  the  angle  with  respect  to  the  vertical  direction.  The  larger  the 
value  of  r,  the  larger  the  anisotropy.  The  other  represents  pancake-like 
distributions: 

(p+l)U 

MD(cos  oc)  =  -*- - sinP  ri.  (11) 

4  up  !! 


Using  Figure  14  of  Dermer  and  Ramaty  (1986),  which  is  for  a  power-law 
electron  energy  spectrum  with  a  power  index  3.5,  I  plot  in  Figure  3  the 
center  to  limb  variations  for  a  downward  beam  distribution  of  r- 3  and  a 
pancake-like  distribution  of  p- 6.  Two  straight  lines  with  a  slope  0.0147  (the 
same  as  that  of  the  straight  line  in  Fig.  2)  are  plotted  for  comparison.  We 
find  that  these  results  are  compatible  with  the  data  shown  in  Figure  2. 
Although  the  shapes  of  the  curves  for  r=3  and  p=6  are  different,  the  data 
points  of  Figure  2  have  too  large  scatter  to  determine  which  case  is  better. 


III.  DISCUSSION  AND  CONCLUSION 


By  using  HXRBS  total  counts  as  calibration,  I  have  studied  the 
directionality  of  y  rays  above  300  keV.  The  normalized  y-ray  limb  brightness 
increases  as  the  heliocentric  angle  9  of  a  flare  increases.  Flares  at  the  limb 
(0=90°)  are  brighter  by  13  times  on  the  average  than  similar  flares  at  the 
central  meridian  (9=15°).  This  behavior  can  be  explained  by  either  pancake¬ 
like  distributions  or  downward  beam  distributions  of  electrons  used  by 
Dermer  and  Ramaty  (1986). 
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There  are  several  ways  to  study  the  directionality  of  energetic 
electrons  accelerated  by  solar  flares.  The  first  is  to  compare  the  number  of 
y-ray  flares  near  the  disk  center  with  that  near  the  limb  (Vestrand  et  al. 
1987;  Rieger  et  al.  1983).  The  second  is  to  study  how  the  spectral  indices  of 
hard  X-rays  and  y-rays  change  with  6  (Bogovalov  et  al.  1935;  Vestrand  et  al. 

1987) .  The  third  is  to  observe  flares  at  different  viewing  angles  (Kane  et  al. 

1988) .  Vestrand  et  al.  have  shown  that  there  are  more  limb  flares  than 
expected  from  an  isotropic  y-ray  emission  pattern.  Rieger  et  al.  have  shown 
that  flares  with  detectable  y-ray  continuum  above  10  MeV  are  all  located 
near  the  limb  (0>6O°).  This  requires  that  relativistic  electrons  with  energies 
>10  MeV  be  distributed  anisotropically  (Dermer  and  Ramaty  1986). 
Bogovalov  ct  al.  have  shown  that  the  spectral  index  of  hard  X-rays  above  50 
keV  decreases  systematically  with  increasing  9.  Vestrand  ct  al.  have  also 
shown  a  similar  behavior  for  the  spectral  index  of  y-rays  above  300  keV. 
These  results  of  different  investigators  using  different  methods  all  point  to 
anisotropic  distributions  of  energetic  electrons. 

In  order  to  deduce  the  distribution  function  of  energetic  electrons 
from  the  number  distribution  of  flares,  one  has  to  perform  detailed 
calculations  of  bremsstrahlung  and  analyze  the  size  distribution  of  flares 
(Petrosian  1985;  Dermer  and  Ramaty  1986;  Dermer  1987).  The  result  of 
Dermer  and  Ramaty  (1986)  shows  that  center-to-limb  variations  larger  than 
10  are  required  to  be  consistent  with  the  excess  y-ray  flares  near  the  limb 
observed  by  Vestrand  et  al.  Petrosian's  (1985)  calculation  also  shows  that 
the  center-to-limb  variation  is  larger  than  10  at  300  keV.  Dermer  and 
Ramaty  have  found  that  the  observed  spectral  hardening  of  limb  flares  is 
consistent  with  r=5  for  a  beam  distribution  and  p>l  for  a  pancake 
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distribution.  However,  stereoscopic  observations  of  Kane  et  al.  show  that  the 
center-to-limb  variation  of  >300  keV  y-rays  is  only  about  2.5  or  less. 

In  the  present  paper,  I  normalized  GRS  fluences  with  HXRBS  total 
counts  to  study  the  anisotropy  of  y-rays.  This  method  is  superior  to  counting 
limb  flares  and  disk  flares  in  the  following  two  scores.  First,  in  this  method 
we  do  not  have  to  consider  the  size  distribution  to  find  out  the  center-to-limb 
variation.  Second,  in  this  method  we  utilize  a  pair  of  values  for  each  flare 
while  in  the  number  counting  method  we  reduce  all  the  data  to  just  two 
numbers.  In  these  respects,  a  study  of  spectral  index  as  a  function  of  0  is 
superior  to  the  number  counting  method. 

The  present  study  shows  that  the  center-to-limb  variation  of  y-rays 
above  300  keV  is  larger  than  10,  which  is  consistent  with  the  results  of 
Petrosian  (1985)  and  Dermer  and  Ramaty  (1986).  Among  recent  studies  of 
solar  y-ray  anisotropy,  only  the  result  of  Kane  et  al.  is  consistent  with  a 
nearly  isotropic  distribution  of  energetic  electrons.  At  present,  the  reason 
for  this  discrepancy  is  not  known.  However,  considering  that  independent 
groups  of  authors  using  different  methods  obtain  consistent  results,  except 
for  Kane  et  al.  (1988),  I  favor  anisotropic  electron  distributions  that  can 
produce  limb  brightening. 

If  we  acquire  y-ray  data  with  larger  dynamic  ranges,  perhaps  with 
Solar  A.  the  planned  Japanese  satellite,  we  will  be  able  to  determine  the 
limb  brightening  factor  with  better  confidence  and  better  accuracy. 

This  research  was  funded  by  ONR  Contract  N00014-85-K-0111,  by 
NASA  Grant  NGL  05-020-272,  and  as  part  of  Solar-A  collaboration  under 
NASA  Contract  NAS8-37334  with  Lockheed  Palo  Alto  Research 


Laboratories.  The  author  thanks  Dr.  C.  Dermer  for  providing  numerical 
results  which  were  published  as  Figures  of  Dermer  and  Ramaty  (1986). 
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FIGURE  CAPTIONS 

Fig.  1 — Correlation  diagram  between  HXRBS  total  counts  and  GRS 
fluences.  The  solid  line  is  the  least  square  fit  to  the  solid  dots  representing 
disk  flares;  the  dashed  line,  "+"  marks  representing  limb  flares.  See  the 
text  for  details. 

Fig.  2--Normalized  GRS  fluence  vs.  heliocentric  angle.  GRS  fluences 
are  normalized  with  respect  to  HXRBS  total  counts.  The  solid  line  is  the  best 
fit  to  the  data. 

Fig.  3— Center-to-limb  variation.  Two  cases  of  Dermer  and  Ramaty's 
calculations  are  shown  together  with  two  straight  lines  having  the  same 
slope  as  the  regression  in  Figure  2, 

Fig.  4— Correlation  diagram  between  HXRBS  total  counts  and  GRS 
fluences  with  the  center-to-limb  variation  removed. 
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